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Composite  high-pressure  hydrogen  storage  vessel  has  been  increasingly  applied  to  the  hydrogen  fuel 
cell  car.  The  design  of  a  composite  vessel  involves  various  integrated  parameters  such  as  the  progressive 
failure  properties,  the  burst  pressure  and  fatigue  lifetime.  The  favorable  combination  of  high  reliability 
and  practicability  of  the  composite  vessel  is  a  challenging  task  from  the  beginning  of  design.  This  paper 
gives  a  comprehensive  review  on  recent  development  of  numerical  simulation  and  optimization  for  the 
designed  composite  vessel.  First,  methods  on  damage  modeling  for  predicting  the  failure  properties  and 
degradation  mechanisms  of  the  composite  vessel  are  reviewed.  Second,  research  on  predicting  the  burst 
pressure  and  lifetime  of  the  composite  vessel  is  reviewed.  The  academic  work  on  the  damage  modeling, 
progressive  failure  analysis  and  finite  element  implementation  which  explains  the  failure  properties  and 
stiffness  degradation  mechanisms  of  the  composite  vessel  is  summarized.  Computational  methods  on  the 
burst  pressure,  the  strength  reliability  and  lifetime  of  the  composite  vessel  are  also  evaluated.  Finally,  ideal 
design  which  aims  to  lessen  the  weight  of  a  composite  vessel  to  the  maximum  extent  under  strength 
and  stiffness  constraints  is  commented.  The  optimization  efficiency  using  different  algorithms  is  also 
comparatively  studied.  The  numerical  simulation  and  optimization  as  important  fundamental  research 
constitute  a  design  platform  for  the  composite  vessel.  It  deserves  pointing  out  the  lightweight  design 
conception  as  a  remarkable  tendency  that  combines  advanced  numerical  methods  and  manufacturing 
technique  develops  rapidly,  commits  to  improving  the  reliability  and  practicability  of  the  composite 
vessel.  It  is  expected  the  lightweight  design  technique  plays  an  increasingly  important  role  in  developing 
the  composite  vessel  as  their  value  is  further  highlighted. 
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1.  Introduction 

The  energy  crisis  and  environment  pollution  all  around  the 
world  enable  it  necessary  to  explore  advanced  energy.  Hydrogen 
energy  is  commonly  considered  as  important  renewable  and  sus¬ 
tainable  energy  in  the  new  century  because  of  a  series  of  advantages 
such  as  cleanness,  high-efficiency  and  fruitful  source  [1-4].  Hydro¬ 
gen  can  be  generated  by  solar  energy,  wind  energy,  biomass  and 
water  [5-7].  Yet,  the  use  of  hydrogen  energy  deals  with  many  tech¬ 
nique  problems  such  as  the  hydrogen  generation,  the  hydrogen 
storage  and  transport,  and  the  hydrogen  combustion.  To  promote 
the  hydrogen  energy  industry  and  potential  hydrogen  economy,  it 
is  critical  to  develop  hydrogen  energy  technique  [8-10]. 

The  hydrogen  fuel  cell  car  is  typically  representative  of  devel¬ 
oping  hydrogen  industry.  Safe,  high-efficiency  and  economical 
hydrogen  storage  technique  is  a  key  to  ensure  favorable  run  of 
hydrogen  energy  car.  Among  many  hydrogen  storage  patterns 
including  high-pressure  gaseous  storage,  cryogenic  liquid  stor¬ 
age  and  chemical  hydrogen  storage,  high-pressure  gaseous  storage 
has  become  the  most  popular  technique  [11-14].  For  example, 
the  hydrogen  fuel  cell  cars  used  in  both  2008  Beijing  Olympic 
Games  and  2010  Shanghai  World  Exposition  used  composite  high- 
pressure  storage  vessel. 

Because  of  the  limited  space  available  for  the  car,  the  basic 
design  requirements  for  the  storage  vessel  are  safety,  reliability 
and  economy.  However,  the  composite  pressure  vessel  may  work 
under  the  high-pressure  and  high-temperature  environment.  This 
not  only  presents  a  strong  challenge  to  the  physical  and  mechan¬ 
ical  performance,  but  also  to  the  reliable  and  economical  design 
about  how  to  achieve  a  perfect  combination  of  safety  performance 
and  low  cost  [15].  Conventional  movable  metal  pressure  vessel  can 
no  longer  be  competent  for  the  rigorous  need  of  high  strength  - 
and  stiffness  -  weight  ratios.  In  this  case,  the  composite  filament 
wound  technology  is  introduced  to  improve  the  storage  efficiency 
[16-18].  To  reach  the  aim  of  large-capacity  hydrogen  storage,  the 
composite  vessel  includes  an  inner  aluminum  layer  and  se  outer 
carbon  fiber/epoxy  composite  wound  layers  in  China  [19].  In  gen¬ 
eral,  the  composites  are  used  to  make  the  pressure  vessel  by  placing 
them  in  different  orientations  for  different  layers  and  in  a  common 
orientation  within  a  layer.  These  layers  are  stacked  in  such  a  way 
to  achieve  high  stiffness  and  strength  [20].  The  designable  proper¬ 
ties  of  the  composite  vessel  enable  it  flexibly  applicable  to  different 
working  conditions  and  environments  by  adapting  the  fiber  wound 
patterns,  the  dome  shapes  and  the  layer  thickness. 

The  design  of  the  composite  vessel  as  a  fundamental  research 
work  relates  the  physical  and  mechanical  properties  of  materials 
to  the  geometry.  Essentially,  two  main  works  are  required  for  the 
design:  one  is  the  prediction  of  burst  pressure  and  another  is  the 
optimization.  The  burst  pressure  denotes  the  limit  load-carrying 
ability  of  the  composite  vessel.  Since  the  failure  mechanisms  of 
composites  are  complicated  from  the  view  of  composite  microme¬ 
chanics,  including  the  matrix  cracking,  fiber/matrix  debonding  and 
fiber  breakage  as  well  as  their  interactions  [21,22],  accurate  pre¬ 
diction  of  the  burst  pressure  is  a  challenging  work.  Besides,  proper 
explanation  of  the  stiffness  degradation  as  the  physical  response  to 
the  continuous  failure  also  needs  special  theory. 

With  the  development  of  computer,  the  finite  element  analysis 
has  become  a  powerful  tool  to  deal  with  the  numerical  problems 
of  complex  structures  [23].  By  combining  the  continuum  damage 
mechanics  (CDM)  and  the  finite  element  analysis  (FEA),  the  pro¬ 
gressive  failure  analysis  of  the  composite  vessel  is  used  to  account 


for  the  stiffness  degradation  and  to  predict  the  burst  pressure 
[24,25].  Especially,  the  FE  technique  that  implements  the  damage 
model  arouses  wide  academic  attention. 

In  terms  of  the  optimization  of  composite  vessel,  several  intelli¬ 
gent  algorithms  that  simulate  the  evolutionary  principles  in  the 
disciplines  of  computational  biology,  physics  and  immunology, 
have  been  used.  The  popular  algorithms  are  the  genetic  algorithm, 
simulated  annealing  and  artificial  immune  system  [26-28].  The 
optimization  attempts  to  improve  the  weight,  strength,  reliability 
and  lifetime  of  the  composite  vessel  by  designing  the  wound  angles 
and  thickness. 

As  the  design  method  and  manufacturing  technique  advance, 
the  composite  vessel  takes  on  the  approach  to  the  lightweight 
development.  The  weight  decrease  has  become  a  technique  bot¬ 
tleneck  which  restricts  the  economy  and  practicability  of  the 
composite  vessel.  In  this  case,  it  is  urgent  that  a  set  of  material- 
structure-processing  integrated  design  and  calculation  method  of 
the  composite  vessel  be  established.  It  is  true  there  are  some  other 
manufacturing  and  test  reasons  affecting  the  performance  of  com¬ 
posite  vessel,  so  what  a  large  role  the  design  plays  may  require 
evaluated.  However,  the  importance  of  numerical  simulation  and 
optimization  is  expected  to  be  further  highlighted  by  the  present 
review. 

2.  Geometry  structure  of  composite  hydrogen  storage 
vessel 

Fig.  1  shows  a  composite  vessel  used  in  the  hydrogen  fuel  car. 
Fig.  2  shows  the  basic  structure  and  wound  pattern  of  composite 
vessel.  The  vessel  includes  an  inner  aluminum  layer  and  several 
outer  composite  layers.  The  aluminum  layer  prevents  the  gas  leak¬ 
age  and  provides  a  mould  for  filament  wound,  and  the  composite 
layers  are  responsible  for  resisting  the  internal  pressure.  For  the 
cylinder  part,  spiral  wound  and  hoop  wound  patterns  are  used.  Yet, 
only  spiral  wound  pattern  appears  at  the  head  part. 

From  the  manufacturing  principle,  the  wound  trace  for  compos¬ 
ite  layers  is  a  combination  of  the  rotation  of  aluminium  liner  and 
the  axial  movement  of  wound  machine.  The  spiral  wound  angle  a0 
at  the  cylinder  is  determined  by  [29] 


Fig.  1.  Lightweight  composite  high-pressure  hydrogen  storage  vessel  used  in  the 
hydrogen  fuel-cell  vehicle. 
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Fig.  2.  Basic  structure  and  wound  pattern  of  composite  vessel:  (a)  hoop  wound  and 
(b)  sprial  wound. 


where  r  is  the  radius  of  the  polar  axis  and  R0  denotes  the  inner 
radius  of  the  cylinder. 

The  spiral  wound  at  the  head  obeys  the  geodesic  path  algorithm, 
marking  the  shortest  distance  between  any  two  points  on  the  head. 
The  spiral  wound  angle  a  at  the  head  changes  from  90°  at  the  polar 
axis  to  oto  at  the  cylinder  [29] 

a  =  arcsin 

The  thickness  H  of  the  composite  layer  at  any  radius  R  of  the 
head  is  calculated  by  [29] 


H  =  h 


(Rp-r2) 
( R 2  -  r2) 


where  h  is  the  wound  thickness  at  the  cylinder.  From  Eqs.  (l)-(3), 
the  polar  radius  r  affects  the  wound  angle  a  and  the  thickness  H  at 
the  head. 


3.  Finite  element  modeling  of  composite  hydrogen  storage 
vessel 

Accurate  and  fast  modeling  of  the  composite  vessel  is  an  impor¬ 
tant  work  in  the  design,  which  involves  many  parameters  such  as 
the  wound  angle  and  thickness.  Parametric  modeling  by  exacting 
these  parameters  as  the  design  variables  facilitates  the  strength 
prediction  and  optimization  for  the  composite  vessels  with  differ¬ 
ent  sizes. 

3.2.  Design  theories 

In  general,  the  composite  vessel  is  considered  as  the  composite 
laminated  structure.  Now,  there  are  two  main  design  theories:  grid 
theory  and  composite  laminate  theory. 

3.2.2.  Grid  theory 

The  basic  assumptions  according  to  the  grid  theory  are:  (1 )  only 
the  longitudinal  carbon  fiber  bears  the  pressure  and  (2)  effects  of 


wound  patterns  are  neglected.  The  grid  theory  can  be  used  to  cal¬ 
culate  the  thickness  of  each  composite  layer  and  the  longitudinal 
in  situ  fiber  strength.  The  strength  value  is  inversely  solved  by  the 
burst  pressure  of  composite  vessel.  Chen  [30]  gave  a  set  of  design 
method  of  the  composite  vessel  based  on  the  grid  theory,  and  the 
total  composite  wound  thickness  hf  is  expressed  as 


/ 


ftp  Pb 

2 Gfjj  COS2  OtQ  ’ 


hf  =  hfa  +  hf0 


3  ftp  Pb 

20fb 


ftp  Pb 

2cfb 


(2  -  tan2  oto ) 


where  hfa  and  hf0  are  longitudinal  and  hoop  wound  thickness,  Pb  is 
the  burst  pressure  of  composite  vessel,  and  cr^  is  the  in  situ  strength 
of  carbon  fiber. 

It  should  be  pointed  out  the  grid  theory  is  merely  an  ideal  design 
method  which  largely  depends  on  the  processing  parameters  and 
the  sizes  of  the  experimental  composite  vessel.  This  shows  seeking 
a  more  efficient  design  theory  is  necessary. 


3. 2 .2.  Composite  laminate  theory 

The  basic  assumptions  of  the  classical  composite  laminate  the¬ 
ory  (CCLT)  are:  ( 1 )  perfect  bonding  appears  at  the  interface  between 
each  layer,  (2)  the  mechanical  properties  of  composite  laminates 
are  substituted  by  those  of  a  middle  plane,  and  (3)  the  normal 
stress  is  neglected  at  the  section  parallel  to  the  middle  plane.  Lif- 
shitz  et  al.  [31  ]  proposed  a  method  using  the  CCLT  to  calculate  the 
stress  and  strain  in  non-symmetric  wound  pressure  vessel  with 
thick  metal  liners.  However,  the  last  two  assumptions  above  do 
not  strictly  hold  for  3D  composite  vessel.  In  this  case,  Parnas  et  al. 
[32]  derived  the  elastic  stress  and  displacement  solutions  for  3D 
cylindrical  composite  laminates  based  on  the  CCLT  and  general¬ 
ized  plane  strain  assumption.  As  the  liner  is  assumed  isotropic 
and  the  composite  layers  are  considered  transversely  isotropic, 
Chapelle  et  al.  [33]  further  derived  the  analytical  solutions  of 
the  stress,  strain  and  displacement  for  3D  cylindrical  composite 
laminates. 

Fig.  3(a)  shows  the  cylindrical  part  of  composite  vessel  and 
Fig.  3(b)  describes  a  representative  volume  element  taken  from 


Aluminum 

liner 


Composite 
layers 


Fig.  3.  (a)  Composite  hydrogen  storage  vessel  (transverse  section)  and  (b)  schematic 
illustration  of  the  relationship  between  the  on-axis  coordinate  (1,2,3)  and  the  off- 
axis  coordinate  ( r,6,z ). 
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Fig.  3(a),  showing  the  principal  direction  ( 1 ,2,3)  of  a  composite  layer 
under  the  cylindrical  coordinate  (r,0,z).  If  the  axial  strains  at  all  lay¬ 
ers  are  assumed  constant  and  the  shear  strains  are  independent  of 
z,  the  strain-displacement  relationship  is  given  by  [33] 


(k)  _  durk)  „(k)  _  urfc)  „(k) 
r  ~  dr  ’  8  ~  r  ’  z 


dUZ}  _  C  v(* 

dz  ~  °’  yz# 


du ^ 
dz 


Yor 

(5) 


where  sz  is  the  axial  strain,  yz0  is  the  shear  strain  and  y0  is  twist 
per  unit  length. 

The  off-axis  stress-strain  relationship  under  the  hygro-thermo- 
mechanics  multiphysics  field  is  expressed  as 
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0's 

_  As  _ 

AC 

(6) 


where  the  symbols  Cy(i,  j  =  1, 2,  3,  6)  are  off-axis  elastic  con¬ 
stants  of  composites.  The  symbols  <rz,  a0 ,  or  and  rz0  are  the  off-axis 
axial,  hoop,  radial  and  shear  stresses,  respectively.  sz,  Sq,  sr  and  yZQ 
are  the  corresponding  strains,  respectively.  AT  and  AC  are  tem¬ 
perature  and  moisture  concentrations,  respectively.  cq(i  =  1, 2, 3,  6) 
and  f3i{i  =  1,  2,  3,  6)  are  the  corresponding  thermal  and  moisture 
expansion  coefficients. 

For  the  kth  layer,  the  equilibrium  equation  under  the  cylindrical 
coordinate  is  given  by  [33] 


r  dr  r 


After  introducing  the  displacement  continuity  conditions,  the 
radial  stress  continuity  conditions  as  well  as  the  axial  equilibrium 
and  zero  torsion  conditions  between  two  neighbouring  layers,  the 
analytical  solutions  for  the  strain,  stress  and  displacement  can  be 
derived  by  Eqs.  (5)-(7). 


3.2.  Finite  element  modeling 


The  finite  element  analysis  and  CCLT  are  combined  to  consti¬ 
tute  the  design  platform  since  the  analytical  solutions  cannot  be 


competent  for  progressive  failure  analysis.  Yan  et  al.  [34]  discussed 
the  finite  element  technique  of  the  composite  vessel  including  finite 
element  modeling  and  commercial  finite  element  codes.  Besides, 
the  reliability,  efficiency  and  cost  as  well  as  the  further  prospect 
are  discussed.  Liu  and  Zheng  [35,36]  established  parametric  finite 
element  model  for  the  cylinder  part  using  ANSYS  finite  element 
code.  Wang  et  al.  [37]  established  the  finite  element  model  of 
the  composite  vessel  with  variable  thickness  and  angles  at  the 
head.  Li  et  al.  [38]  established  the  design  model  of  the  compos¬ 
ite  vessel  using  ANSYS-APDL  (Ansys  Parametric  Design  Language) 
and  VC++  language.  However,  the  design  platform  has  not  been 
completely  established  since  ANSYS  software  does  not  provide  an 
efficient  modeling  module  for  composites  until  now.  Based  on  the 
composite  micromechanics  and  modified  CCLT,  Antunes  et  al.  [39] 
established  the  finite  element  model  of  the  composite  vessel  using 
ABAQUS  and  Algor  finite  element  codes.  Park  et  al.  [40]  calculated 
the  filament  wound  patterns  using  semi-geodesic  path  equations 
for  an  arbitrary  surface,  and  established  the  finite  element  model  of 
the  composite  vessel  using  ABAQUS  finite  element  code.  Further, 
Kim  et  al.  [41]  used  the  semi-geodesic  path  algorithm  to  calcu¬ 
late  possible  wound  patterns  by  considering  the  windability  and 
slippage  between  the  fiber  and  the  mandrel  surface.  Yet,  these 
work  concentrates  on  the  specific  composite  vessel  and  does  not 
perform  parametric  modeling.  Recently,  we  established  the  para¬ 
metric  finite  element  model  by  combining  the  ANSYS  and  ABAQUS 
finite  element  codes.  Fig.  4  shows  the  flow  chart  of  finite  ele¬ 
ment  modeling  of  composite  vessel,  and  Fig.  5  shows  a  90°  mesh 
model.  A  subroutine  is  written  using  ANSYS-APDL  to  establish  the 
geometry  and  mesh  model,  where  some  characteristic  parame¬ 
ters  such  as  the  composite  lay-up  thickness,  the  cylinder  inner 
radius  and  thickness,  the  liner  thickness  as  well  as  the  mouse 
radius  and  thickness  are  extracted  as  the  design  variables.  Then, 
the  mesh  model  is  imported  into  ABAQUS  by  the  ANSYS/ABAQUS 
file  transformation  interface.  After  that,  a  local  cylindrical  coor¬ 
dinate  system  (r,0,z)  is  defined  in  ABAQUS  and  the  composite 
lay-up  module  is  employed  to  assign  the  lay-up  angles  and  mate¬ 
rial  properties  for  composite  elements.  Fig.  6  shows  schematic 
stress  distributions  under  working  pressure  for  a  composite  ves¬ 
sel.  The  wound  angles  at  the  cylinder  are  90°,  -90°,  a0,  -Qfo,  90°, 
-90°,  — g?o,  90°  and  -90°  in  turn.  It  can  be  seen  that  the  maxi¬ 
mum  principal  stress  appears  at  the  joint  of  the  head  and  cylinder 


Characteristic 

parameters 


Wound 

parameters 


Layup  angle 

VLayupV 

thickness 

Cylinder 

length 


Finite  element 
analysis 
using 
ABAQUS 

t 


Geometry 

Parameters 


Cylinder 
inner  radius 

Liner 

thickness 

Mouse 
inner  radius 


Finite 

element 

modeling 

using 

ANSYS- 

APDL 


Mouse  length 


Fig.  4.  Flow  chart  of  finite  element  modeling  of  composite  vessel. 
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Fig.  5.  Parametric  finite  element  model  of  composite  vessel. 


for  the  spiral  wound  layer,  in  good  agreement  with  Yang’s  results 
[29]. 

4.  Failure  analysis  and  strength  prediction  of  composite 
hydrogen  storage  vessel 


According  to  Lemaitre  and  Chaboche  [51,52],  Ristinmaa  and 
Ottosen  [53],  Schapery  [54],  Murakami  and  Kamiya  [55],  Hayakawa 
et  al.  [56]  and  Simo  et  al.  [57],  the  damage  constitutive  modeling 
based  on  the  CDM  is  summarized.  For  the  damaged  materials,  the 
Cauchy  stress  tensor  a  can  be  substituted  by  the  nominal  stress 
tensor  a 


a  =  (1  —  D)a 


where  0<D<1.  D  =  0  represents  the  perfect  materials  and  D  =  1 
denotes  the  completely  damaged  materials. 

The  Helmholtz  free  energy  per  unit  mass  xfr  for  elastic  materials 
under  isothermal  conditions  is  written  as 


pf  =  pfe(eij-eI!j,Dij)  +  pfd(K )  (10) 

where  xf/e  and  x//d  are  the  free  energy  which  represents  the  elastic 
deformations  and  damage  hardening,  respectively,  k  is  an  internal 
variable  and  p  is  the  density. 

The  thermodynamic  conjugate  forces  (Yy,  B)  corresponding  to 
the  internal  variables  (Dy,  k)  are  expressed  as 


G  ij  ~  P 


dx/f 

W 


p 


dx/f 

W 


(ii) 


The  Clausius-Duhem  dissipation  inequality  according  to  the 
CDM  is  written  as 


y = -py + au  £ij  >  o 


(12) 


Currently,  finite  element  technique  has  been  widely  used  to 
predict  the  failure  properties  and  strengths  of  composites.  In  the 
following,  some  recent  research  on  the  failure  analysis  and  strength 
prediction  of  the  composite  vessel  is  reviewed. 


where  y  is  the  power  of  dissipation  due  to  damage.  If  the  dissipation 
y  reaches  the  maximum,  the  damage  evolution  law  is  given  by 


8Fd 

w 


k  =  -X 


ddF d_ 
dB 


(13) 


4.1.  Damage  constitutive  modeling 

With  the  accelerating  application  to  the  hydrogen  fuel  cell  car  in 
China,  the  composite  vessel  may  be  subjected  to  the  hygro-thermo- 
mechanics  environment  arising  from  the  subtropical  climate.  The 
aging  of  epoxy  matrix  leads  to  the  stiffness  and  strength  degrada¬ 
tion,  thereby  affecting  the  lifetime  of  composite  vessel.  In  this  case, 
it  is  significant  to  explore  the  failure  properties  using  progressive 
failure  analysis. 

Strictly  speaking,  there  are  two  distinct  failure  mechanisms 
for  the  composite  vessel:  intralaminar  damage  and  interlaminar 
delamination.  From  the  burst  experiments  of  composite  vessels, 
the  intralaminar  damage  instead  of  interlaminar  delamination 
affects  the  mechanical  performance  of  the  composite  vessel.  This 
is  because  two  neighbouring  composite  layers  still  keep  normally 
pressing  and  tangentially  slipping  despite  a  little  interlaminar  fail¬ 
ure  appearing  before  the  burst  [42].  Thus,  only  the  intralaminar 
failure  is  generally  considered. 

Since  Kachanov  [43]  proposed  the  CDM  to  study  the  creep  rup¬ 
ture  of  metals,  progressive  failure  analysis  using  the  CDM  has 
been  widely  used  to  predict  the  stiffness  degradation  and  damage 
evolution  properties  of  composites  by  defining  phenomenological 
damage  tensors  and  damage  evolution  laws  [44-50].  A  schematic 
3D  damage  model  from  Eq.  (6)  under  the  hygro-thermo-mechanics 
field  is  expressed  as 


After  introducing  three  damage  potential  functions  Fd(Yy,  D,  B) 
for  three  intralaminar  failure  mechanisms  for  composites,  the  con¬ 
sistency  parameter  Xd,  the  damage  variable  D  and  elastic  matrix  Q 
are  obtained  as  the  functions  of  stress  and  strain. 

4.2.  Failure  criteria  for  composites 

Failure  shows  the  degradated  ability  to  resist  the  external 
load.  The  failure  properties  of  a  composite  vessel  depend  on  the 
geometry  sizes,  the  material  features,  the  loading  history  and  envi¬ 
ronment.  Over  the  past  five  decades,  the  popular  failure  criteria 
are  the  maximum  stress,  Hashin,  Hoffman,  Hinton,  Christensen, 
Rotem,  McCartney,  Yamada-Sun,  Chang-Chang,  Tsai-Hill,  Tsai-Wu, 
Puck  and  Huang’s  criteria  [58-70].  In  terms  of  the  tensor  notation, 
the  generalized  Tsai-Wu  criterion  is  written  as 

FjCTj  +  Ffj Gj Gj  +  Fij^GfGjGj^  >  1 ,  i,j,  k  =  1,2,  . . . ,  6  (14) 

where  F  is  strength  coefficient  tensor. 

Such  as  the  Tsai-Wu,  Tsai-Hill  and  Hoffman  failure  criteria 
identify  the  failure  of  composites,  but  they  cannot  identify  the  fail¬ 
ure  modes.  In  contrast,  the  maximum  stress,  Hashin,  Christensen, 
Rotem,  McCartney,  Yamada-Sun,  Hinton,  Chang-Chang,  Puck  and 
Huang  criteria  can  identify  the  failure  modes  by  defining  the  fail¬ 
ure  criteria  for  the  tensile  and  compressive  failure  as  well  as  the 
fiber  breakage,  matrix  cracking  and  delamination. 
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where  dfi=f  m,  s)  are  the  damage  variables  representing  the  Hinton  et  al.  [60]  explored  existing  failure  theories  for  corn- 

fiber  breakage,  matrix  cracking  and  fiber/matrix  interface  failure,  posites  by  launching  a  world-wide  failure  exercise,  and  tested 
respectively.  some  leading  failure  theories  by  comparing  experimental  data. 
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Fig.  6.  Schematic  stress  distributions  at  (a)  liner,  (b)  hoop  wound  layer  and  (c)  spiral  wound  layer  for  a  designed  composite  vessel. 


Nahas  [71],  Tsai  [72]  and  Sleight  [73]  performed  extensive  liter¬ 
ature  surveys  on  the  failure  theories  of  composites  and  classified 
the  failure  criteria  into  two  types:  the  non-interactive  failure 
criteria  and  interactive  failure  criteria.  The  non-interactive  fail¬ 
ure  criterion  is  defined  as  one  without  interactions  between  the 
stress  or  strain  components.  These  criteria  compare  the  individ¬ 
ual  stress  or  strain  components  with  the  corresponding  material 
allowable  strengths,  and  indicate  the  type  of  failure  modes.  The 
maximum  stress,  Hashin,  Christensen,  Rotem,  McCartney,  Yamada- 
Sun,  Chang-Chang,  Puck  and  Huang’s  criteria  belong  to  this  type. 
In  contrast,  the  interactive  failure  criteria  consider  interactions 
between  stress  components.  These  theories  use  a  polynomial  based 
on  the  material  strengths  to  describe  a  failure  surface.  The  Tsai-Wu, 
Tsai-Hill  and  Hoffman  failure  criteria  belong  to  this  type.  Currently, 
there  is  also  an  evolving  trend  to  develop  the  fracture  energy-based 
failure  criteria  to  predict  the  failure  of  composites  using  such  as 
the  cohesive  fracture  energy-based  criterion  [74,75]  and  finite  ele¬ 
ment  failure-based  criterion  [76,77].  However,  accurate  prediction 
of  fracture  properties  of  composites  requires  advanced  experimen¬ 
tal  technique. 

4.3.  Progressive  failure  analysis 

Essentially,  progressive  failure  analysis  of  the  composite  ves¬ 
sel  is  to  implement  the  proposed  damage  constitutive  model, 
the  failure  criteria  and  damage  evolution  law.  Garnich  and  Akula 
[78],  Orifici  et  al.  [79],  and  Liu  and  Zheng  [24]  gave  reviews 
on  the  current  methods  on  the  progressive  failure  analysis  and 
finite  element  implementation  of  composites.  The  finite  ele¬ 
ment  equation  under  the  hygro-thermo-mechanics  field  is  given 
by 

(15) 


where  I<=  f  yBTCdBdV  is  the  stiffness  matrix  of  a  finite  element, 

Fs  =  fsNTFsdS,FT  =  f  vBTCda  A TdV and  Fc  =  f  vBTCdjS  AC  dV are  ele¬ 
ment  forces  arising  from  the  surface  load,  thermal  expansion  and 
moisture  expansion,  respectively.  B  =  LNT  is  strain  matrix,  L  is  dif¬ 
ferential  operator  and  N  is  shape  function. 

The  flow  chart  for  the  progressive  failure  analysis  of  the  com¬ 
posite  vessel  using  finite  element  method  is  illustrated  in  Fig.  7,  as 
summarized  by  five  points  [73]:  (1)  for  each  load  step,  the  finite 
element  analysis  is  performed  and  the  on-axis  stresses/strains  at 
each  element  are  obtained,  (2)  the  stresses/strains  at  each  ele¬ 
ment  are  compared  with  the  material  allowable  values  and  used 
to  determine  whether  some  elements  fail  according  to  the  failure 
criteria.  If  no  failure  is  detected,  the  applied  load  is  increased  and 
the  analysis  continues,  (3)  if  some  elements  fail,  the  damage  vari¬ 
ables  are  updated,  and  the  element  stiffness  constants  are  degraded 
according  to  the  damage  model.  After  that,  the  equilibrium  of  the 
composite  vessel  is  re-established  using  the  modified  stiffness.  This 
adjustment  accounts  for  the  nonlinear  solution  due  to  the  stiffness 
change,  (4)  the  calculations  are  performed  continuously  under  the 
same  load  until  there  is  no  failed  element  anymore,  (5)  the  iterative 
process  repeats  until  the  catastrophic  failure  showing  the  burst  of 
the  composite  vessel. 

Chang  [80]  and  Ju  et  al.  [81]  performed  the  initial  and  progres¬ 
sive  failure  analysis  of  the  composite  vessel  based  on  the  CCLT. 
However,  the  analysis  is  only  limited  to  the  analytical  solution. 
Further,  Perreux  and  Thiebaud  [82],  Liu  and  Zheng  [35]  per¬ 
formed  the  progressive  failure  analysis  of  the  composite  vessel 
by  proposing  damage  models  and  developed  the  corresponding 
finite  element  technique.  The  numerical  convergence  problem  is 
solved  by  introducing  the  arc-length  algorithm  or  viscous  sta¬ 
bilization  method  [83,84].  The  fracture  localization  problem  is 
addressed  by  adding  a  characteristic  parameter  into  the  damage 
model  to  eliminate  the  mesh  effect  [85].  However,  to  develop  the 


Ku  =  Fs  +  Fj  +  Fc 
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Fig.  7.  Flow  chart  of  progressive  failure  analysis  of  composite  vessel. 


finite  element  technique  for  the  failure  analysis  under  the  hygro- 
thermo-mechanics  multiphysics  field  may  still  require  a  long 
way. 

Since  the  failure  of  composite  is  the  macroscopic  representation 
of  the  microscopic  fiber/matrix  failure  evolvement,  it  is  insufficient 
to  explore  the  failure  properties  of  composite  vessels  merely  by 
using  the  CDM  theory.  In  the  past  two  decades,  there  were  continu¬ 
ous  efforts  for  developing  hierarchical  and  multiscale  finite  element 
methods  to  perform  the  progressive  failure  analysis  of  composites. 
The  multiscale  damage  model  at  the  microscopic  and  macroscopic 
scales  for  the  composite  vessel  is  shown  in  Fig.  7.  Square  and  hexag¬ 
onal  cell  models  are  established  at  the  microscopic  scale  based  on 
the  observation  of  carbon  fiber/epoxy  composites  using  scanning 
electron  microscope.  Thus,  the  macroscopic  damage  tensors  above 
can  be  represented  by  microscopic  damage  information.  He  and 
Wang  [86]  studied  the  transverse  matrix  cracking  and  stress  con¬ 
centration  for  the  composite  vessel  by  multiscale  failure  analysis. 
Haj-Ali  and  Muliana  [87]  proposed  an  integrated  micromechan¬ 
ics  framework  for  the  nonlinear  analysis  of  composites.  However, 
these  work  is  only  limited  to  the  simple  composite  structures.  Fur¬ 
ther,  Blassiau  et  al.  [88]  established  the  microscopic  cell  models 
on  the  fiber/matrix  load  transfer  and  the  multiscale  model  on  the 
progressive  failure  analysis  of  the  composite  vessel.  However,  the 
failure  properties  of  the  microscopic  epoxy  matrix  was  not  con¬ 
sidered  and  the  damage  initiation  information  for  the  fiber/matrix 
interface  cannot  be  acquired.  With  the  development  of  computer, 
the  computational  ability  for  the  multiscale  analysis  linking  the 
microscopic  to  macroscopic  failure  analysis  is  largely  improved. 
Potential  difficulty  lies  how  to  reveal  all  microscopic  failure  mech¬ 
anisms  in  a  realistic  way  which  may  be  elusive  of  engineering  tests 
[74]. 


4.4.  Burst  pressure  prediction 

Burst  pressure  of  the  composite  vessel  is  the  ultimate  aim  of 
progressive  failure  analysis.  Mao  et  al.  [89]  estimated  the  burst 
pressure  based  on  the  statistical  approach  and  experiments.  Cohen 
et  al.  [90,91]  investigated  the  effects  of  manufacturing  and  design 
variables  such  as  the  material  parameters  and  wound  parameters 
on  the  strength  and  stiffness  of  the  composite  vessel  by  exper¬ 
iments.  Hwang  et  al.  [92,93]  performed  the  probabilistic  failure 
analysis  and  experimental  tests  to  predict  the  size  effect  on  the 
burst  pressure  of  composite  vessels  using  Monte  Carlo  method. 
However,  no  information  on  the  progressive  failure  properties  of 
composites  is  obtained  from  these  work.  Sun  et  al.  [94]  performed 
the  nonlinear  finite  element  analysis  to  predict  the  burst  pressure 
of  the  composite  vessel  using  the  maximum  stress  failure  criterion 
and  a  stiffness  degradation  model.  Onder  et  al.  [95]  explored  the 
influence  of  temperature  and  wound  angle  on  the  burst  pressure  of 
the  composite  vessel  by  the  modified  CCLT  and  experiments.  How¬ 
ever,  they  defined  only  the  constant  degradated  stiffness  regardless 
of  the  evolvement  of  damage  variables.  Based  on  the  modified  CCLT 
and  CDM,  Liu  and  Zheng  [35,96]  proposed  damage  models  and  finite 
element  algorithm,  respectively,  to  predict  the  burst  pressure  of  the 
composite  vessel  by  considering  the  progressive  failure  properties 
and  nonlinear  damage  evolvement.  Prediction  of  the  burst  pres¬ 
sure  is  essentially  related  to  the  numerical  convergence  problem 
in  the  finite  element  analysis.  According  to  the  design  requirement 
of  the  composite  vessel,  only  little  matrix  cracking  under  work¬ 
ing  pressure  is  permitted.  As  the  pressure  increases,  progressive 
interface  debonding  and  matrix  cracking  appear.  Only  when  a  large 
number  of  fiber  breakage  suddenly  appears  is  the  composite  vessel 
considered  to  burst. 
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5.  Optimal  design  of  composite  hydrogen  storage  vessel 

The  aim  of  design  for  the  composite  vessel  is  to  reach  the 
best  match  between  the  reliability  and  practicability.  Advanced 
design  technique  can  shorten  the  duration  and  save  the  cost.  Now, 
the  ISO  (International  Standardization  Organization)  and  European 
Union  have  issued  the  standards  for  the  design  of  composite  ves¬ 
sel  [97,98].  Unluckily,  there  is  still  no  mature  national  standard  of 
composite  vessel  until  now  in  China  since  most  of  the  design  and 
manufacturing  depend  on  the  experience  and  experiment  [99,19]. 
In  this  case,  developing  a  set  of  efficient  design  method  is  nec¬ 
essary.  Now,  two  main  design  methods  for  the  composite  vessel 
sprout:  safety  coefficient-based  design,  and  reliability  and  lifetime- 
based  design,  which  are  considered  as  the  improvement  over  the 
conventional  experiential  design. 

5.2.  Safety  coefficient-based  design 

According  to  Leung  [100],  the  working  pressure  of  composite 
vessel  is  taken  as  the  burst  pressure  divided  by  the  safety  coeffi¬ 
cient.  However,  the  difficulty  lies  in  how  to  select  a  proper  safety 
coefficient  which  yet  depends  on  the  design  and  manufacturing 
levels.  The  CGH2R  Draft  Revision  10  of  European  Union  and  the 
ISO/DIS  15869.2  specify  the  safety  coefficients  of  2.35  and  2.25, 
respectively  [97,98].  Then,  the  optimal  design  attempts  to  reach 
the  minimum  weight  and  highest  strength  under  the  constraint  of 
safety  coefficient.  Tummala  et  al.  [101]  performed  the  multiobject 
optimization  of  the  composite  vessel  on  the  pressure,  temperature 
and  weight.  Krikanov  [102]  proposed  a  new  optimization  method 
to  design  the  composite  vessel  with  the  minimum  mass  under 
strain  and  strength  constraints.  Especially,  Liang  et  al.  [103]  car¬ 
ried  out  the  optimal  design  of  dome  contours  for  the  composite 
vessel  subjected  to  geometric  limits  and  wound  condition.  Wang 
and  Zheng  [37]  proposed  a  mathematical  design  model  based  on 
the  finite  element  analysis.  However,  the  optimization  efficiency 
above  is  low  since  only  simple  optimization  methods  are  used.  In 
this  case,  searching  advanced  intellectual  algorithms  becomes  nec¬ 
essary.  Ghiasi  et  al.  [104]  compared  the  optimization  efficiency  of 
various  methods  such  as  the  gradient-based  method,  direct  search 
method,  genetic  algorithm  and  simulated  annealing.  Messager  et  al. 
[105]  developed  a  genetic  algorithm-based  optimization  method 
to  maximize  the  stability  buckling  pressure  of  the  composite  ves¬ 
sel.  Kim  et  al.  [106]  proposed  an  optimal  design  method  of  the 
composite  vessel  by  associating  the  finite  element  analysis  with 
the  genetic  algorithm.  Liu  et  al.  [107,108]  proposed  an  adaptive 
genetic  algorithm  and  an  immune  algorithm,  respectively,  to  per¬ 
form  the  optimal  design  of  the  composite  vessel  which  aims  to 
reach  the  minimum  weight  under  the  constraints  of  safety  coef¬ 
ficient  and  processing  parameters.  The  parametric  model  for  the 
composite  vessel  is  shown  in  Fig.  8.  The  intellectual  optimization 
algorithms  show  higher  convergence  speed  and  precision  than  the 
gradient-based  method. 

5.2.  Reliability  and  lifetime-based  design 

For  a  designed  composite  vessel,  economy  and  practicability 
are  based  on  safety  and  reliability.  The  reliability  and  lifetime  of 
the  composite  vessel  relate  to  many  factors  such  as  the  material 
parameters,  wound  patterns  and  manufacturing  defects.  It  is  not 
always  feasible  to  expect  to  use  a  single  safety  coefficient  to  pro¬ 
vide  the  guarantee  for  any  occasional  accident.  Thus,  it  is  necessary 
to  develop  the  reliability  and  lifetime-based  methods  to  perform 
the  probabilistic  design.  Compared  with  the  deterministic  design 
method,  the  probabilistic  design  method  may  avoid  costly  over- 
design  and  conservative  manufacturing. 
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Fig.  8.  Optimization  model  of  composite  vessel. 

Now,  the  reliability  and  lifetime-based  design  conception  has 
been  successfully  applied  to  the  composite  vessel.  Beakou  et  al. 
[109,110]  optimized  the  wound  angle  of  cylindrical  composite 
laminates  using  reliability  methods  for  the  composite  structure. 
Gomes  et  al.  [Ill]  performed  the  optimization  of  laminated  com¬ 
posites  under  the  reliability  constraint  using  genetic  algorithm  and 
neural  network.  However,  the  geometry  model  in  there  work  is 
much  simpler  than  the  composite  vessel.  Further,  Shen  et  al.  [112] 
launched  the  reliability-based  probabilistic  design  of  the  compos¬ 
ite  vessel  by  combining  the  statistical  theory  and  experiments. 
However,  the  probabilistic  design  is  only  limited  to  the  analytical 
solution  using  the  grid  theory.  Cai  et  al.  [113]  developed  a  set 
of  reliability-based  design  method  of  the  composite  vessel  and 
performed  the  sensitivity  analysis  on  the  influence  of  statistical 
variables  on  the  safety  factors  under  external  pressure.  Based  on 
the  progressive  failure  analysis,  Liu  and  Zheng  [114]  performed 
the  strength  reliability-based  design  of  the  composite  vessel  using 
Monte  Carlo  simulation  and  response  surface  method,  in  which 
the  lay-up  thickness  and  angles  of  composite  layers  are  consid¬ 
ered  as  the  stochastic  variables.  Besides,  Wang  [115]  compared 
several  existing  lifetime  prediction  technique  of  the  composite 
vessel  under  static  and  cyclic  internal  pressure.  Camara  et  al. 
[116]  studied  the  strength  reliability,  the  failure  probabilities  and 
the  lifetime  of  the  composite  vessel  by  multiscale  modeling  and 
composite  micromechanics.  Bie  et  al.  [117]  proposed  a  theoretical 
model  to  predict  the  fatigue  lifetime  of  the  composite  vessel  by 
combining  the  micromechanics  and  CDM.  However,  the  research 
above  did  not  touch  upon  the  lifetime-based  design.  Comond  et  al. 
[118]  proposed  a  set  of  fatigue  lifetime-based  design  method  of 
the  composite  vessel  based  on  the  modified  CCLT  and  the  CDM. 
However,  the  design  for  the  head  part  was  not  considered  in 
their  work.  Further,  Bertin  et  al.  [119]  explored  the  effects  of 
wound  patterns  on  the  fatigue  lifetime  of  the  composite  vessel 
under  thermomechanical  cyclic  loading  by  experiments.  It  can  be 
concluded  the  strength  reliability  and  lifetime-based  material- 
structure-processing  integrated  design  technique  will  become  an 
important  tendency  for  the  design  of  composite  vessel.  The  best 
design  method  is  to  set  up  the  high-fidelity  numerical  simulation 
and  advanced  optimization  platform.  But,  the  design  must  also 
consider  too  many  factors  including  the  statistical  properties  of 
fiber  strengths,  complex  geometry  structure,  and  efficiency  and 
stability  of  numerical  calculations.  For  the  probabilistic  design,  the 
prediction  of  the  burst  pressure  and  fatigue  lifetime  of  composite 
vessel  is  main  limitation  which  otherwise  needs  more  funda¬ 
mental  research.  After  all,  the  complicated  failure  mechanisms  of 
composites  are  still  beyond  our  understanding  until  now  [120]. 
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6.  Further  prospect  of  composite  hydrogen  storage  vessel 

In  recent  years,  the  composite  pressure  vessels  develop  rapidly 
in  the  field  of  hydrogen  storage.  Academic  and  engineering  atten¬ 
tions  are  paid  to  the  design  and  manufacturing  technique  of  them. 
With  the  development  of  computer-aided  design  technique,  the 
composite  vessel  offers  designers  much  design  space.  The  design  is 
to  reach  the  ideal  match  among  the  reliability  and  lifetime,  weight 
and  cost. 

It  is  commonly  recognized  the  composite  vessel  should  take  on 
the  lightweight  road  [121].  The  lightweight  design,  manufacturing, 
test  and  evaluation  technique  guarantee  the  safe  and  economical 
application  of  composite  vessel.  The  ideal  design  method  is  to  set  up 
a  set  of  material-structure-processing  integrated  design  platform, 
and  to  solve  the  numerical  simulation  and  optimization  problems 
efficiently.  This  is  not  easily  realized  in  view  of  the  complicated 
relationships  between  the  design  requirements  and  design  param¬ 
eters. 

The  composite  vessel  is  an  important  application  of  such  as 
the  advanced  carbon  fiber/epoxy  composites  to  the  area  of  hydro¬ 
gen  energy.  The  complicated  failure  mechanisms  and  degradation 
mechanisms  are  distinct  characteristic  of  composites  although  they 
exhibit  high  stiffness-  and  strength-density  ratios.  Further,  the 
hygro-thermo-mechanics  environment  is  also  a  large  problem  in 
the  design  of  the  composite  vessel.  More  fundamental  research  on 
the  failure  properties  of  composites  under  the  multiphysics  field 
remains  to  be  performed. 
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